Superconductivity in MgB 2 and TaB 2 : A Full-Potential Electronic Structure 

Comparison 
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We present a comparison of electronic structure of MgBi and Ta£?2, the two new superconduc- 
tors, as well as VB2, calculated using full-potential, density- functional-based methods in P6/mmm 
crystal structure. Our results, described in terms of (i) density of states (DOS), (ii) band-structure, 
and (iii) the DOS and the electronic charge density in a small energy window around the Fermi en- 
ergy Ef- In particular, the charge density around Ef in MgP>2 and TaB^ show striking similarity 
as far as the B plane is concerned. A comparison of their band-structures, coupled with /-character 
analysis, indicates that TaBi has substantially more p-character than VB2 along A — L and H — A 
directions near Ef- 
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With the discovery of superconductivity in MgB 2 pi , 
many Boron compounds such as BeB2, LiB^, ATB2 B] 
and transition-metal borides |?|-[| with ^M_B 2 - S tructure 
are being reexamined for possible superconductivity. The 
search for new MgB 2 -\ike superconductors is motivated 
not only by the desire to discover new superconductors 
but also with the hope that it will be helpful in unrav- 
elling the nature of interaction responsible for supercon- 
ductivity in these materials. In this context, the recently 
reported superconductivity in TaB 2 Q with a supercon- 
ducting transition temperature T c ps 9.5K, if confirmed 
may lead to a better understanding of supercon- 
ductivity in MgB 2 and MgB 2 -like superconductors such 
as TaB 2 . 

The continuing experimental |^-|| and theoretical 
P-|l3| efforts at characterizing and understanding the na- 
ture of superconductivity in MgB 2 have made substan- 
tial progress. However, a complete understanding of su- 
perconductivity in MgB%, and for that matter in TaB 2 , 
is still lacking. It is interesting to note that the supercon- 
ducting transition temperatures in MgB 2 and TaB 2 are 
39if and 9.5K, respectively, while VB 2 has been found 
not to superconduct down to 0A2K Thus a the- 

oretical study of those aspects of electronic structure of 
MgB 2 , TaB 2 and VB 2 , which could be responsible for 
superconductivity or lack thereof in these materials, may 
provide some clues as to why TaB 2 is superconducting 
but VB 2 is not. This, in turn, may improve our under- 
standing of superconductivity in MgB 2 itself. 

In this paper we present a full-potential electronic 
structure study of TaB 2 and VB 2 in PQ/mmm crys- 
tal structure using density-functional-based methods. 
For comparison with MgB 2 , we use our fully-relaxed, 
full-potential electronic structure results for MgB 2 in 
PQ/mmm crystal structure fl3|| . We analyze our re- 
sults in terms of (i) density of states (DOS), (ii) band- 
structure along symmetry directions, and (iii) the DOS 
and the electronic charge density in a small energy win- 
dow around the Fermi energy, Ep, in these systems. 



We have used our own full-potential program as well as 
LMTART package jl4| to calculate, self-consistently, the 
electronic structure of TaB 2 and VB 2 using the experi- 
mental lattice constants. These calculations have been 
used to study the band-structure along symmetry di- 
rections and site- and symmetry-decomposed densities 
of states. For studying the charge density in a small 
energy window around Ep we have used the Stuttgart 
TB LMTO package |5). The fully-relaxed, full-potential 
electronic structure of MgB 2 was calculated as described 
in Ref. @ 

Based on our calculations, described below, we find 
that in TaB 2l within a small energy window around Ep, 
the charge density in the B plane is still two-dimensional 
and it is associated with the p x ,y orbitals as is the case 
for MgB 2 |h],[Ii],[13| . The corresponding charge density 
in VB 2 does not have this character. A comparison of 
band-structures of MgB 2l TaB 2 and VB 2 along symme- 
try directions reveals that the point is well below Ep 
in both TaB 2 and VB 2 . However, based on our under- 
standing of superconductivity in MgB 2 |p|-|l3|, we find 
that the 7-th band, specially around K and A points 
with substantially more p-character in TaB 2 than in VB 2 
, may play an important role in making TaB 2 a super- 
conductor. 

Before describing our results in detail, we provide some 
of the computational details of our calculations. The 
charge self-consistent full-potential LMTO calculations 
were carried out with the generalized gradient approxi- 
mation for exchange-correlation of Perdew et al. [lq] and 
484 fc-points in the irreducible wedge of the Brillouin 
zone. The basis set used consisted of s, p, d and / orbitals 
at the Ta (V) site and s, p and d orbitals at the B site. 
The potential and the wavefunctions were expanded up 
to Imax — 6. The input to the tight-binding LMTO cal- 
culations, carried to charge self-consistency, were similar 
to that of the full-potential calculations except for us- 
ing spherically symmetric potential and the space-filling 
atomic spheres. 
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In Table I we give the lattice constants used for calcu- 
lating the electronic structure of MgB2, TaB 2 and VB 2 . 
Note that the lattice constants for TaB 2 and VB2 are 
taken from experiments 0] , while that for MgB 2 are ob- 
tained theoretically as described in Ref. |U5). In Table I 
we also show the calculated plasma frequencies and the 
total density of states per spin, n(Ep), at Ef for these 
compounds. 





a (a.u.) 


c (a.u.) 


< (eV) 


(eV) 


n(E F ) 


MgB 2 


5.76 (5.834) 


6.59 (6.657) 


7.04 


6.77 


4.70 


TaB 2 


5.826 


1.0522 


9.01 


9.71 


6.57 


VB 2 


5.665 


1.0196 


5.62 


6.60 


9.24 



We show in Fig. 1 the total density of states for MgB 2 , 
TaB 2 and VB 2 calculated using the full-potential LMTO 
method at the lattice constants given in Table I. The 
gross features of the DOS of TaB 2 and VB 2 are simi- 
lar if one takes into account the differences arising out 
of the change in the d-bandwidth as one goes from V 
to Ta. The comparison with MgB 2 must be done keep- 
ing in mind the differences in the atomic electronic con- 
figurations between Mg (3s 2 ) and Ta (5d 3 6s 2 ) and V 
(3d 3 is 2 ). As a result the e?-DOS at Ef changes from 0.9 
st/Ry at the Mg site to 8.8 st/Ry at the Ta site, and 
finally to 16.8 st/Ry at the V site in the respective di- 
borides. More importantly, we hnd that the p-DOS at 
the B site in MgB 2 . TaB 2 and VB 2 decreases steadily 
from 3.36 st/Ry in MgB 2 to 1.55 st/Ry in TaB 2 and 
to 0.58 st/Ry in VB 2 . The decrease in the availability 
of p-electrons around Ep may result in the decrease of 
T c in TaB 2 and no superconductivity in VB 2 . The DOS 
of TaB 2 compares reasonably well with the cluster cal- 
culations of [0, however, the location of Ep is clearly 
different. 

In the conventional, BCS-type superconductors the 
DOS within an interval of ±hwp> (u>D=Debye frequency) 
at Ep is crucial for superconductivity. Since the B par- 
tial DOS is more relevant for superconductivity in the 
diborides, we have plotted in Fig. 2 the symmetry- 
decomposed DOS at the B site within a small energy 
interval around Ep in MgB 2 , TaB 2 and VB 2 . From 
Fig. 2 we find that in TaB 2 the DOS at Ep with x (y) 
symmetry is half of the corresponding DOS in MgB 2 , 
while in VB 2 it is more than a factor of four less than in 
MgB 2 . Since the electrons in the B plane are expected to 
play a crucial role in determining superconducting prop- 
erties, it is not surprising to find significant differences 
in the DOS with x (y) symmetry in MgB 2 , TaB 2 and 
VB 2 . 

To study the differences in the band-structures of 
MgB 2 and the transition metal diborides as well as be- 
tween TaB 2 and VB 2 , in Fig. 3 we have plotted the 
band-structures along the symmetry directions for TaB 2 
and VB 2 . The presence of d-electrons in the transition 
metal diborides implies more d-like bands below Ep in 
comparison to MgB 2 . As a result, the point is well 
below Ep in both TaB 2 and VB 2 . However, the 7-th 



band at K point, with significant p z character, in TaB 2 
is below Ep while it is above Ep in VB 2 . Our /-character 
analysis also shows that the 7-th band, which crosses Ep 
along A — L and H — A directions, has substantial p x and 
p z characters in TaB 2 but only a very small p character 
in VB 2 . Based on our understanding of superconductiv- 
ity in MgB 2 |||l(| , it seems reasonable to state that the 
7-th band, specially around K and A points, may play 
an important role in making TaB 2 a superconductor. 

To further illustrate the differences between MgB 2 , 
TaB 2 and VB 2 , we have calculated their electronic 
charge densities within a 5 mRy energy window around 
Ep using the TB LMTO method with the lattice con- 
stants as given in Table I. In these calculations the atomic 
sphere radii were adjusted to minimize the discontinuity 
in the Hartree potential across the atomic spheres which 
have been found to give reliable results ]is| . In Fig. 4, 
we show the charge density calculated within a 5 mRy 
energy window around Ep for TaB 2 and VB 2 . By com- 
paring the charge density of TaB 2 with that of MgB 2 as 
given in Ref. p3| , we find that to some extent the two- 
dimensional character of the charge density due to p x>y 
orbitals in the B plane remains in tact in TaB 2 . In VB 2 , 
however, the decrease in the lattice constants a and c 
has completely destroyed the two-dimensional nature of 
the charge density in the B plane, These differences may 
account for no superconductivity in VB 2 , 

In conclusion, we have presented the results of our 
full-potential electronic structure calculations for MgB 2 , 
TaB 2 and VB 2 in PQ/mmm crystal structure using 
density-functional-based methods. We have analyzed our 
results in terms of the density of states, band-structure 
along symmetry directions, and the DOS and the elec- 
tronic charge density in a small energy window around 
the Fermi energy in these systems. In TaB 2 , we find that 
within a small energy window around Ep, the charge den- 
sity in the B plane is still two-dimensional and it is asso- 
ciated with the Px.y orbitals as is the case for MgB 2 . The 
corresponding charge density in VB 2 does not have this 
character. By comparing their band-structures we find 
that the 7-th band, specially around K and A points with 
substantially more p-character in TaB 2 than in VB 2l 
may play an important role in making TaB 2 a super- 
conductor. 

The author would like to thank A. V. Mahajan for 
helpful discussions. 
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FIG. 2. The partial B density of states around Fermi en- 
ergy calculated at the lattice constants as given in Table I 
using the full-potential LMTO method. 
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FIG. 4. The isosurfaces of charge density within a 5 mRy 
energy window around the Fermi energy in the primitive cell 
oiTaB2 (left) and VB2 (right), calculated at the lattice con- 
stants as given in Table I using the TB LMTO method. The 
values of the isosurfaces (8,6,4,2,1) decrease as one moves 
away from the B plane. The unit for charge density in the 
volume is 10~ 4 [e/(a.w.) 3 ] and the view is off-diagonal at an 
angle of 30°. 
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